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ABSTRACT 

In this paper we develop the new semi-analytical model of a tidally perturbed or tidally 
disrupted star proposed recently by two of us. This model is effectively a one dimensional 
Lagrangian model and it can be evolved numerically much faster than the conventional 
3D models. A self-consistent derivation of dynamical equations of the model is performed 
and several important theorems about the dynamics of the model are proved without any 
particular assumption about the equation of state of the stellar gas. The dynamical equations 
are solved numerically for the case of an n = 1.5 polytropic star evolving in the relativistic 
tidal field of a 10'' M© Kerr black hole. Some results of these calculations are compared with 
the results of calculations based on finite-difference 3D simulations. The comparison shows 
a very good agreement between both approaches to the problem. Then we show that the 
strength of the tidal encounter depends significantly on the relative orientation of the orbital 
angular momentum of the star and the spin of the black hole. 
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2 P. B. Ivanov et al. 
1 INTRODUCTION 

This paper continues studying a new dynamical model of a star evolving under the influence 
of a tidal field. This semi-analytical model has been proposed by two of us in a recent paper 
(Ivanov & Novikov, 2001, hereafter IN). It allows researchers to calculate the outcome of a 
strong tidal encounter of a star with a source of a tidal field much faster than the standard 
three-dimensional approach, and evolves the star under the influence of a tidal field during 
a much longer time. On the other hand, testing of the model for the case of a Newtonian 
tidal field of a point mass has shown that the new model gives very good quantitative 
agreement with the results of 3D simulations. Therefore, the new model could be used in 
a situation where the formulation of a problem demands many different computations of 
the tidal encounter events with different parameters, or for calculating the stellar evolution 
in a complicated tidal field, and where the present day 3D simulations cannot be used due 
to problems with computational time or other problems A natural example of such a 
problem is an attempt of surveying the parameter space of the problem of tidal disruption 
of stars by a Kerr black hole. 

The main feature of the new model consists of an assumption about the motion of the 
stellar gas perturbed by the action of a tidal field. Namely, it is assumed that different 
mass shells of the star always keep the shape of ellipsoids during the evolution of the star 
in the tidal field. [J] This assumption allows us to reduce the complicated non-linear three 
dimensional dynamics of the stellar gas to an effectively one-dimensional Lagrangian nu- 
merical scheme. The dynamical equations of our model are derived from the so-called virial 
relations written for each mass shell (see the next Section), and form a set of non-linear 
one-dimensional partial differential equations of hyperbolic type coupled with the tidal field. 
They depend on time and a Lagrangian variable which could be represented by the mass 
enclosed within a particular shell or the radius of the shell in the unperturbed spherical state 
of a star. 

There is one fundamental drawback in the variant of the model studied by IN. Namely, 
the 'thermal terms' (i.e. the terms determined by the gas energy and pressure) in the virial 

* See IN for an overview of works on tidal encounters and astrophysical applications. 

t Therefore, the model is a direct generalisation of the so-called afHne model of the star (Carter & Luminet, 1982, 1983, 1985, 
Lattimer & Schramm, 1976) which has been intensively used for investigation of the tidal encounters before the development of 
modern 3D computational methods. However, contrary to the affine model, the different elliptical mass shells evolve differently, 
with different parameters and orientations. 
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relations were treated by IN in an approximate manner. This led to unphysical behaviour 
of the model, the mass shells corresponding to different Lagrangian radii were allowed to 
intersect each other during the evolution of the star. Therefore, the variant of IN was not 
suitable for the study of the density, pressure and velocity distribution within the star. 
However, the distribution of these quantities represents a significant interest in the problems 
connected with the problem of tidal disruption. For example, one would need to know these 
quantities for a study of the subsequent evolution of the gas liberated from the star after a 
fly-by of a black hole and the formation of an accretion disc (or torus). The main purpose of 
this paper is to resolve this difficulty of the model. In this paper, we calculate the 'thermal 
terms' exactly and derive the dynamical equations of our model in a self-consistent way (see 
the next Section). Then, we apply our advanced variant of the model to the problem of a 
fiy-by of a n = 1.5 polytropic star around a Kerr black hole. We test the model against three 
dimensional calculations made by Diener et al (Diener et al, 1997) for the same problem 
and the same parameters. We find very good agreement between our calculations and the 
calculations based on the 3D approach. 

Our paper is organised as follows. We derive the dynamical equations of our model in 
the next Section. In Section 3 we discuss the results of numerical simulations. Discussion 
and conclusions are presented in Section 4. The paper is written in a self-consistent way, 
and all important relations are derived in the text without referring to IN. 

Following IN we use an unusual summation convention: summation is performed over 
all indices appearing in our expressions more than once but summation is not performed 
if indices are enclosed in brackets. Bold letters represent matrices in abstract form. All 
indices can be raised or lowered with help of the Kronecker delta symbol, but we distinguish 
between the upper and lower indices in order to enumerate the rows and columns of matrices, 
respectively. 

2 DERIVATION AND ANALYSIS OF THE DYNAMICAL EQUATIONS 

We derive and analyse dynamical equations of our model in a manner similar to that was 
described in our previous paper (IN). However, as we mention in the Introduction, we do 
not use simplifying assumptions concerning "thermal terms" in our equations (for the exact 
definition of the thermal terms, see equations (15-19) below). Also, we prove several impor- 
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tant theorems about the dynamics of our model without any particular assumptions about 
the equation of state of the stellar gas. 

At first, we would like to introduce coordinate systems and several useful kinematical 
quantities. We use two different coordinate systems: a) Cartesian coordinates associated 
with a locally inertial frame centered at the star's geometrical center (we call those below "the 
Eulerian coordinates" of the gas element); b) Cartesian coordinates Xq of the gas elements 
in an unperturbed spherical state of the star (say, before the star is deformed by the tidal 
field) . By definition, these coordinates are not changing during the evolution of any particular 
gas element, and therefore we call them below "the Lagrangian coordinates". As we have 
mentioned in the Introduction, we assume that the star consists of elliptical shells, and these 
shells are not deformed during the evolution of the star. This assumption allows us to write 
the law of transformation between the Lagrangian and Eulerian coordinates in the form: 



where ro = yxoiXQ is the Lagrangian radius of a particular shell, and Cq = x^/ro are 
direction cosines in the Lagrangian space (eojeg = 1). We represent the position matrix T 
and its inverse S as a product of two rotational matrices A and E, and a diagonal matrix 
B: 

Ti = A]BlEJ^ = a,A]E^, = a^'A>El, 2 
where = a(i)6^, and ai are the principal axes of the elliptical shell. The Jacobian 



D = |f4| of the mapping between the Lagrangian and Eulerian spaces can be written as 



D{Xo) = - 



'0 

where 

g = |T| = 0x0203, 4 

is the determinant of the matrix T, and the symmetric matrix i?'" determines a local shear 
and a change of volume of the neighbouring shells: 

^'" = ^(sL(Tr)' + 5;^(Tno, 5 

with the primes standing for differentiation with respect to ro 1^. The matrix R can be 

t Note a useful formula for the averaged value of the Jacobian D (IN): D = J dClD = Here integration is performed 

over a unit sphere in the Lagrangian space. 
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represented in terms of its eigenvalues fm and the rotational matrix O: 

-R'" = fmO^jD"^. 6 

As it follows from the law of mass conservation, the evolution of the gas density p is deter- 
mined by the evolution of D\ 

where po('^o) is the gas density in the unperturbed spherical state of the star. If one of 
the eigenvalues fm goes to zero, the density p(t, x*) goes to infinity at a certain value of 
Xq. This physically corresponds to the intersection of two neighbouring shells. However, we 
assume that pressure forces can always prevent the shells from intersecting, and therefore 
the eigenvalues fm are always positive. 

Similar to IN, we start the derivation of the dynamical equation of our model from 
the integral consequences of the exact hydrodynamical equations: the equation of energy 
conservation and the so-called virial relations. We write the energy conservation equation in 
adiabatic approximation, thus neglecting the energy transfer between neighbouring shells, 
the entropy generation due to nuclear reactions, viscosity, etc., 

J d^x{pv^/2 + e)+V} = - J dS,{pv') + j d^xipCijv'x^), 8 

Here f* is the velocity of the gas element, v = y/viV^, p is the pressure and e is the energy 
density per unit volume. V stands for the potential energy of the star. The matrix Cj 
represents the tidal tensor, and therefore it is symmetric and traceless. The virial relations 
have the form: 

d^x{px^v') = J c/^x(pfV) + 5^' j d^xp 
" dSiix'^p) +V'''+ f d^x{pCix''x^), 9 



where V^^ is the so-called potential energy tensor: 



T /" /" f ly'^ rr''^\ ( ly'^ ly'^ 



V"' = — / d'x / d'xip(x')p(x])^ , i^. 10 

2J J \x — Xi\ 

Obviously, we have V = V^^. 

Now we substitute the evolution law (1) in equations (9) and calculate all terms in these 
equations. Then, differentiating the result with respect to the Lagrangian mass coordinate 

pro 

M{ro) = 47r / po{ri)dri 11 



we obtain the dynamical equations of our model (we use below the mass M(ro) of the gas 
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inside the shell of radius tq as a new Lagrangian coordinate instead of tq). Analogously we 
obtain a differential form of the law of energy conservation from the integral form (8). 

We calculate explicitly "gravitational" parts of equations (8) and (9) (the potential en- 
ergy and potential energy tensor) making an additional simplifying approximation. Namely, 
we assume that the gravitational force acting on the gas near the shell with some Lagrangian 
radius ro is equivalent to the gravitational force of a uniform density ellipsoid with a mass 
equal to the part of the star's mass within that shell. The principal axes of that ellipsoid 
coincide with the principal axes of the shell, and the density is averaged over the volume 
enclosed in the shell. Under this assumption the "averaged" potential energy tensor V^^ has 
the form 

= -- [ dMGMMA''^^^, 12 
2 J 9 

and the "averaged" potential energy V = V^^ is 

V= — / dMGM^^. 13 
27 g 

The dimensionless quantities Dj have been described in e. g. Chandrasekhar, 1969, p. 41. 
They have the form: 

^ r°° du 

Dj = 9 . . 2 I V 1^ 
Jo A(a^ + u) 



where A = ^ (af + u) (al + u) (al + u) . 

For the "thermal" terms 11*^ = 6^^ J d^xp — J dSiX^p in equation (9), we obtain 

ff^ = S''' I dM{^) - AngSiTl'P^^ 15 
J p 



±fdn^-, 16 



where 

^ p' 477 J p 

and 

P^-' = i- I dnpeie^,, 17 

and the integration is performed over a unit sphere in Lagrangian space. Analogously, the 
energy term / d^xe in equation (8) has the form: 

.3. r ....e. 



d-^xe = / rfM(-), 18 
J p 

where (^) = ^ / dfl^, and the surface term / dSi{pv^) has the form: 

J dSiipv') = A-ngP^^S^tl 19 
The calculation of other terms in equations (8), (9) is straightforward. Differentiating equa- 
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The new model of a tidally disrupted star: further development and relativistic calculations 7 
tion (8) with respect to M and taking into account (13), (18,19), we obtain: 

+ 3(-) - 2^'^^— > = -^'^-^dM^aP^'S^ti) + qT/T/. 20 
Equation (20) is a differential form of the law of energy conservation. It is analogous to 
equation (22) of IN. Analogously, differentiating equation (9) and taking into account (12) 
and (15), we obtain: 

T] = 3Si(^ - l2nSi^{gSlT::P^-} - ^jT.a.D.E^^ + CIT^ . 21 

Equations (21) are the dynamical equations of our model. They are analogous to equations 
(23) of IN. 

Equation (20) must follow from equations (21). To prove this, we contract both sides 
of equations (21) with the velocity matrix over all indices and subtract the result from 
equation (20). The remainder is separated into gravitational and thermal parts. As it was 
shown by IN, the gravitational part is reduced to the identity: 

— (^^) + ^ ^ ^ = 0. 22 
dt g g 

For the thermal part we have: 



^(i) + 3if(^) + ^p'" 

dt p p por^ 



+ 3H{q + ^P'^Ri^ = 0, 23 



where the expansion rate is defined as 

H=\slTi = \i^ + '-^ + '-^), 24 
o o ai a2 CI3 

and we change variables in the last term according to rule (11). Now we are going to prove 
that equation (23) follows directly from the first law of thermodynamics written in the 
adiabatic approximation: 

d ,e ^ d ,1^ 
dt p dt p 

For that, we differentiate equations (3), (7) with respect to time to find: 

D = 3HD + ^eoiConi?'", 26 



and 

p D 

Substituting (26), (27) in equation (25), and using (7), we obtain 

^(i) + 3i/^ + ^P.-' 
dt p p por^ 

Integrating (28) over the solid angle Q, we obtain (23). Therefore, equation (23) is equivalent 



27 



+ 3H^ + ^PeyoRin = 0. 28 
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to the zeroth moment of equation (25): 

d ,e . r dVL d ,1. ^ 

-ri-) + / -^V^i-) = 0. 29 
dVp' J A-n^dVp' 

Now we are going to show that the quantities 

XAM) = Tir; - Tifi 30 

are exact integrals of motion of the dynamical system (21). For that, we contract the left 
and right hand sides of equations (21) with and take the antisymmetric part of the result. 
We obtain 

The quantity in braces in (31) can be transformed to 

pnj j^kn _ pnkpjn 22 

To prove that (32) is equal to zero, it is sufficient to note that both symmetric matrices P 
and R can be diagonalized by the same orthogonal transformation (defined with the help of 
the matrix O, see equation (6)). Therefore, the commutator (32) of P and R must be equal 
to zero. As was discussed by IN, the quantities Xjk{M) describe the circulation of the fluid 
along our elliptical shells, and therefore the conservation law 

Xjk{M) = 0. 33 

is a direct consequence of the conservation of circulation in our model. 

Let us discuss the law of conservation of angular momentum. In our model the angular 
momentum tensor L^"^ can be expressed as 

L^' = - f dMiT^tl - T^f^). 34 
3 ^ 

Let the quantity l^^ be the angular momentum tensor density per unit of mass: l^'^ = -^L^'^. 
Then, one can easily obtain the evolution law for l^"^ from equations (21): 

i'' = in^{gP'''{SiT^ - Sir!;)} + l{Ci,T-T} - C^-T^. 35 
The first term in (35) describes the transfer of angular momentum between neighbouring 
shells due to pressure The second term is obviously the tidal torque term. The quantity 
in braces is equal to zero in the center and also at the boundary of the star. Therefore, if 
the tidal torque is absent, the angular momentum is conserved. 

§ Note that transfer of angular momentum due to self-gravity is absent. Obviously, this is related to our simplified description 
of the self-gravity forces in our model. 



© 1994 RAS, MNRAS 000, 0-?? 



The new model of a tidally disrupted star: further development and relativistic calculations 9 
3 NUMERICAL WORK 

For our numerical work we would like to choose a simple polytropic model of the star with 
the polytropic index n = 1.5 and the specific heat ratio 7 = 5/3. The star has the radius 
R^, and mass M^, equal to the radius and mass of the Sun. It is assumed that the star 
is moving along a marginally bound orbit in the field of the Kerr black hole. The same 
problem has been discussed by Diener et al for the case of rather weak tidal interaction, and 
we use the results of this work for testing our model in the relativistic tidal field. We use 
the natural 'stellar' units in our calculations and representation of results: the dimensionless 



time r = \J^-, the radius R = r/R^,, the mass coordinate x = M/M^, energy E = 
and specific angular momentum L = L/^/GM^R^. 

As was pointed out by IN, in the non-relativistic approximation the whole problem can 
be described by only two parameters: the polytropic index n, and the quantity 



M, Rl 



\M,hRr 

where Mbh is the black hole mass and Rp^ is the pericentric separation from the black hole 
calculated in the non-relativistic approximation. This quantity has been introduced by Press 
and Teukolsky in the linear theory of tidal perturbations. The smaller this quantity is, the 
stronger the tidal encounter will be. For the relativistic field the situation is much more 
complicated. For a fixed stellar model, the problem must be parametrised by the ratio of 
the star's mass to the black hole mass in order to specify the importance of relativistic 
corrections. The problem also depends on the dimensionless rotational parameter a of the 
black hole 1^. The marginally bound star's orbit can be specified by its angular momentum 
Lorb and the Carter integral Q. In this paper we would like to consider the most interesting 
case of equatorial orbits, and therefore we set Q = 0. We use the dimensionless orbital 
angular momentum Lorb = Lorb/rgC, where rg = ^^^5^. Instead of using of 77^^, Frolov et al 
(1994) proposed to use a more relevant quantity 

V = yinr{Rp/R;'f'\ 37 

where Rp is the minimal radial distance from the black hole calculated in the relativistic 
approach. Note that the corresponding dimensionless quantity Rp = Rp/ rg can be expressed 
only in terms of Lorb provided the rotational parameter a is specified. 

For our calculations we use an explicit Lagrangian numerical scheme which is similar to 

^ The dimensionless rotational parameter a is determined by the black hole angular momentum J and its mass M: a = Q^^i ■ 
© 1994 RAS, MNRAS 000, 0-?? 
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Figure 1. Evolution of total (solid curves), thermal (long dashed curve), gravitational (short dashed curve) and kinetic (dot- 
dashed curve) energies with time. All energies are measured with respect to their equilibrium values. The dotted curve represents 
the evolution of total energy of the gravitationally bound part of the star in the model of Diener et al. 

what was used by IN. The stabihty criterion of our scheme is discussed in Appendix 1. The 
essential differences of the new numerical scheme from the scheme of IN are discussed in 
Appendix 2. The results of our calculations are presented in Fig 1-9. 

In Figures 1 — 4 we show the time dependence of different energies, the angular momen- 
tum, central density and mass lost by the star. We choose M^h = 1.0853-10^0, a = 0.9999, 
7] = 1.6486 and Lorb = 2.72945. These parameters correspond to the Diener et al model 5 
which has been intensely investigated in 3D finite difference simulations. Therefore, we use 
this calculation for testing our model and numerical scheme. In Fig. 1 we show the time de- 
pendence of the total energy of the whole star (the upper solid curve) and the total energy 
of gravitationally bound debris (the lower solid curve) [jj]. The dot-dashed curve corresponds 
to the kinetic energy of the star, the dashed curve corresponds to the gravitational energy 
of the star and the long-dashed curve corresponds to the thermal energy of the star. In 
general, all curves look very similar to the corresponding curves calculated by IN in the 
non-relativistic approximation. The dotted curve shows the time dependence of the total 
energy of the gravitational bound debris calculated by Diener et al in the 3D calculations. 
One can see from this Figure that this curve almost coincides with our curve. 

In Fig. 2 the total angular momentum of the star (solid curve), the angular momentum of 
the gravitationally bound debris (dashed curve) and the same quantity calculated by Diener 

II The gravitationally bound debris is defined as the part of the star where the sum of kinetic and gravitational energies is less 
than zero. 
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Time 



Figure 2. Evolution of the component of angular momentum perpendicular to the equatorial plane of the black hole. 
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Figure 3. Evolution of the central density (expressed in units of the central density of the unperturbed star) as a function of 
time. 

et al (dot-dashed curve) are shown as a function of time. The total angular momentum of the 
star grows monotonically with time, and is significantly larger than the angular momentum 
of the gravitationally bound debris at the end of the calculations. The angular momentum 
of gravitationally bound debris calculated in the 3D simulations is close to our quantity for 
r < 6. Then, a sharp decrease of the angular momentum is observed and at the end of the 
calculation the angular momentum of Diener et al is significantly smaller than our quantity. 
The reason for this behaviour of the angular momentum found in the 3D simulations is not 
clear to us, and therefore we cannot comment on this deviation. 

In Fig. 3 the time evolution of the the central density (expressed in units of the central 
density of the unperturbed star) is shown. The solid curve corresponds to our model and 
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Figure 4. The amount of mass (expressed in units of the stellar mass) lost by the star during the tidal encounter as a function 
of time. 

the dashed curve corresponds to the 3D simulations. We can see again that these two curves 
are very similar. 

In Fig. 4 we show the amount of mass lost by the star (expressed in units of the stellar 
mass) with time. The solid curve and the dashed curve correspond to our model and the 3D 
simulations, respectively. The asymptotic value of the mass loss is almost the same in both 
cases and is about 0.5. However, we would like to note that there is a significant ambiguity 
in determining the gravitationally bound gas in the 3D simulation (Diener et al). If one 
considers all the gas elements which left the computational domain with velocities less than 
the escape velocity as still being present in the debris, the mass loss would be significantly 
less with an asymptotic value ~ 0.3. In general, we think that our model shows very good 
similarities to the 3D simulations for the parameters chosen in this computation. 

Figure 5 shows the time evolution of the velocity field and density contours in the orbital 
(XOY) plane for the same model. Axis (OX) is directed opposite to the black hole during 
the pericenter passage. Time r = corresponds to the time of the pericenter passage. In 
this figure one can see that the density contours lose their spherical shape with time and 
elongate. The size of the outermost contour at r = 5 is about 10 times larger than at 
r = —2. In the beginning, the contours expand in the direction of the black hole, but as the 
star approaches the pericenter of its orbit, they start to lag behind. The lag of the innermost 
contours is slightly less than that of the outer ones. The distribution of the velocities in the 
star is represented by arrows whose lengths are proportional to the velocity magnitude. 

Now let us discuss some simple properties of the parameter space of the problem. Figures 

© 1994 RAS, MNRAS 000, 0-?? 
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Figure 5. The time evolution of the velocity field and density contours in the orbital X-Y plane for the rj = 1.6486 encounter. 
The spatial scales are expressed in units of the unperturbed stellar radius i?* . The large arrow points toward the black hole. 
The innermost contour corresponds to 95% of the central density. At each subsequent contour density decreases by a factor 
of 1.5 . Upper left-hand panel: t = —2 {R/Rp = 1.5), Upper right-hand panel: r = {R/Rp = 1), Bottom left-hand panel: 
r = 3 {R/Rp = 1.9), Bottom right-hand panel: r = 5 {R/Rp = 2.73). The small arrows show the direction and the relative 
magnitudes of the velocities. 

(6-9) show dependencies of the mass lost by the star, the total energy contained in the 
gravitationally bound debris and the angular momentum of the debris on the orbital angular 
momentum Lorb for two values of the rotational parameter of the black hole: a = 0.9999 
and a = 0. The mass of the black hole and other orbital parameters are the same as in 
the previous calculations. The largest orbital angular momentum corresponds to rj = 2, the 
smallest angular momentum corresponds to the total disruption of the star (we denote the 
respective value of rj as the 'critical' rjcr)- 
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Figure 6. The mass lost by the star after a fly- by of the black hole is shown as a function of the absolute magnitude of the 
dimensionless orbital momentum. The case Mf^f^ = 1.0853 ■ IO^Mq is considered. The solid and dashed curves correspond to 
a = 0.9999, positive and negative orbital angular momenta, respectively. The dotted curve correspond to a = 0. 




1.9 2.0 



Figure 7. Same as Figure 6, but the mass lost by the star is shown as a function of the parameter rj. 

Figures (6,7) show the dependence of the mass lost by the star on the absolute magnitude 
of the orbital angular momentum and parameter 77, respectively. The solid and dashed curve 
correspond to positive and negative angular momenta and are calculated for a = 0.9999. The 
dotted curve is calculated for a = 0. It is seen from the figures that the stars with negative 
orbital angular momenta are disrupted much more effectively than the stars with positive 
angular momenta 0. It is obvious that similar curves calculated for smaller rotational pa- 
rameters of the black hole (but with the same mass) must lie between the solid and dashed 



** Note that this effect has been discussed by Beloborodov et al (1992) in the framework of a rather naive criterion for tidal 
disruption. 
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Figure 8. The energy of gravitationally bound debris is shown as a function of the absolute magnitude of the dimensionless 
orbital angular momentum. The solid and dashed curves correspond to a = 0.9999, positive and negative orbital angular 
momenta, respectively. The dotted curve correspond to a = 0. 



0.30 - 




Orbital angular momentum 

Figure 9. The angular momentum of gravitationally bound debris is shown as a function of the absolute magnitude of the 
dimensionless orbital angular momentum. The solid and dashed curves correspond to a = 0.9999, positive and negative orbital 
angular momenta, respectively. The dotted curve correspond to a = 0. 

curves. We have rjcr = 1-5 for a = and rjcr = 1.28, 1.69 for a = 0.9999 and positive and 
negative orbital angular momenta, respectively. 

In Figure 8 we show the dependence of the total energy of the gravitationally bound 
debris calculated after the fly-by of the black hole on the value of the dimensionless orbital 
angular momentum. Similar to Figures 6,7 representing the mass loss. Figure 8 shows that 
the stars moving on orbits with negative orbital angular momentum are perturbed more 
effectively than the stars moving on orbits with positive angular momentum. It is interesting 
to note that the asymptotic value of the total energy in the limit of small orbital angular 
momentum (corresponding to full disruption of the star) is a nonzero quantity which depend 
neither on the spin of the black hole nor on the sign of the orbital angular momentum (see 
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Figure 8). This could be explained as follows: after the fly-by the stellar gas leaves the star 
with almost zero specific total energy (i. e. parabolic velocities). So, the out-flowing gas 
does almost not carry specific energy, and the energy of the gravitationally bound part of 
the star is conserved. In Figure 9 we show the component of the angular momentum of the 
star perpendicular to the orbital plane after the fly-by as a function of the dimensionless 
orbital angular momentum. Contrary to the total energy, the angular momentum of the 
gravitationally bound debris is a non-monotonic function of the dimensionless orbital angular 
momentum. 

4 DISCUSSION 

In this paper we construct a self-consistent variant of the new model of a tidally perturbed or 
tidally disrupted star proposed by Ivanov and Novikov 2001. The model allows researchers 
to calculate the outcome of the tidal disruption event with the help of a one-dimensional 
Lagrangian numerical scheme. Therefore, it is much faster than the conventional numerical 
3D approach, and it could be evolved for much longer time. We use the model in numerical 
calculations of the tidal interaction of an n = 1.5 polytropic star with a Kerr black hole of 
mass 1O^M0. We compare the results of our calculations with the results of finite difference 
"iD calculations of the same problem and find a very good agreement between them. Then, we 
consider dependencies of the main characteristics of the tidally perturbed star after a fiy-by of 
the black hole in the equatorial plane on the value of the orbital angular momentum. We find 
that the stars with negative orbital angular momentum are perturbed more effectively than 
the stars with positive orbital angular momentum. We also briefly discuss the dependence 
of the outcome of the tidal encounter on the spin of the black hole. 

As it was demonstrated in the present work (see also IN), the model gives results which 
in certain cases coincide almost completely with results of 3-D calculations. On the other 
hand, the dynamical equations of our model cannot be reduced to the exact hydrodynam- 
ical equations. Therefore, a natural question arises: why is the agreement between the two 
approaches so good? The possible explanation might be as follows. The key assumption of 
our model consists in using elliptical shells for the description of the shape of the star evolv- 
ing under the influence of the tidal fleld. It seems that the quadrupole dependence of the 
fleld of tidal forces on the angular coordinates and the special algebraic properties of the 
tidal tensor could justify such assumptions at least for large scale hydrodynamical motions 
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induced in the star. This could answer quahtatively the question why different elements of 
the star which are not in causal contact evolve in such a way that the elliptical form of the 
shells is always maintained. 

Now let us discuss the problems of our model. At first we discovered in our numerical 
calculations our numerical is scheme is slowly unstable for stellar models with a more sharp 
density contrast. For example, in the case of an n = 3 polytrope it takes several characteristic 
'stellar' times for the instability to halt the computations. Since our numerical scheme has 
been written in a rather naive manner, we expect that a more sophisticated numerical scheme 
(e.g. an implicit scheme) could resolve this difficulty. There is a more fundamental problem 
of our model. We managed to obtain the distribution of pressure and density across the star 
in a simple form only for a polytropic star. Therefore, it is not clear for us how to use our 
model for the case of a more realistic stellar gas [j^. 

In the present paper we do not make an attempt to comprehensively survey the parameter 
space of the relativistic tidal problem and to calculate cross-sections of different kind. This 
will be treated in future work. 

ACKNOWLEDGEMENTS 

We are grateful to J. C. B. Papaloizou and A. G. Polnarev for very useful discussions, and 
Martin Goetz for very helpful comments. This work has been supported in part by RFBR 
grant 00-02-16135, in part by the Danish Research Foundation through its establishment of 
the Theoretical Astrophysics Center, and in part by the Danish Natural Science Research 
Council through grant N 9701841. 

REFERENCES 

Beloborodov, A. M., lUarionov, A. F., Ivanov, P. B., Polnarev, A. G., 1992, MNRAS, 259, 209' 
Diener, P., Frolov, V. P., Khokhlov, A. M., Novikov, I. D., Pethick, C. J., 1997, ApJ, 479, 164 
Carter, B., Luminet, J.-P., 1982, Nature, 296, 211 
Carter, B., Luminet, J.-P., 1983, A&A, 121, 917 
Carter, B., Luminet, J.-P., 1985, MNRAS, 212, 23 

Chandrasekhar, S., 'Ellipsoidal Figures of Equilibrium' 1969, (New Haven: Yale Univ. Press) 
Frolov, V. P., Khokhlov, A. M., Novikov, I. D., Pethick, C. J., 1994, ApJ, 432, 680 
Ivanov, P. B., Novikov, I. D., 2001, ApJ, 549, 467 
Lattimer, J. M., Schramm, D, M., 1976, ApJ, 210, 549 



tt Note that the variant of the model considered by IN is free from this difficulty. 
© 1994 RAS, MNRAS 000, [|-?? 



18 P. B. Ivanov et al. 

Press, W. H., Teukolsky, S. A. 1977, ApJ, 213, 183 



APPENDIX A: PROPAGATION OF SMALL PERTURBATIONS AND THE 
TIME STEP CONSTRAINT 

In a standard approach the time step constraint must follow from a stability analysis of 
numerical schemes. However, the standard stability analysis of our numerical scheme is 
rather complicated and therefore we do not use it in our paper. Following IN, we constrain 
our time step by the condition 

5t = . Al 

^smax 

Here a < 1 is a parameter. Csmax is greater than or equal to the velocity of propagation of 
a small perturbation (with respect to the mass coordinate) Cg calculated in analytical linear 
approximation: Csmax > c^. For the stellar gas we assume the equation of state of an ideal 
gas with polytropic index 7. 

To estimate Csmax we decompose our dynamical variable T in a background part and a 
perturbation: T = Tq + t, where the perturbation t is assumed to be of standard oscillatory 
form: 

t = te'^MM)_ ^2 

For the velocity Cg we have 

Ul 

The dynamical equation for the perturbation follows directly from equations (21) and 

(17). Neglecting the dependence of the background quantities on the mass coordinate and 

taking into account only the second term on the right hand side of equation (21), we have: 
n2 n2 q2 

with the symmetric tensor 

J 47r (egeg/??)(^+') 

Substituting (A2) in (A4), one can obtain a set of algebraic equations. Then, the usual 
compatibility condition gives the value of the velocity c^. However, this is too complicated 
for our purposes. To estimate the upper limit of the velocity, we can use an upper limit 
estimate for the tensor H^i^"^-. 

jjljkn ^ _^ (^^Ij^kn ^ ^Ik^ln ^ ^In^lk-^^ ^5 

J- 3 /mm 
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where fmin is the minimal eigenvalue of the matrix R. Substituting (A5) in (A3), and using 
(A2), we have: 

where 



4 = SWj. 



and 



J min 

Let US assume the matrix a to be symmetric: a = a-^. In this case both matrices a^- and 5*' 5*/ 
can be diagonalized by the same transformation. Taking into account that the quantities 
a^ are the eigenvalues of the matrix S\Sl, equation (A6) can be reduced to the form: 

A 

Ci = — a-2(2Q + ^Q), A9 

i=l 

where the quantities q are the eigenvalues of the matrix a. Equation (A9) has nontrivial 
solutions if and only if 

h («?^ - 2) 

where x = c^/A. Obviously, equation (AlO) gives an implicit dispersion relation. A simple 
analysis of equation (AlO) shows that all roots of (AlO) are larger than 

/ = 3(1 + 1 + 1), All 

af ai ai 

and we have 



Therefore, we use the condition 



Cs > JfA. A12 



St = ^ 413 



in order to constrain our time step. 



APPENDIX B: THE NUMERICAL SCHEME 

The numerical scheme used in the computations is very similar to what was used by IN 
and we address the reader to Appendix B of IN for a comprehensive description of the 
scheme. The main difference is determined by the fact that we now calculate the exact 
distributions of pressure and density across the star contrary to the approximate treatment 
of these quantities by IN. These quantities determine the thermal terms (16,17) which in 
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turn determine the action of pressure forces in our model. In this Appendix we outhne the 
numerical evaluation of the thermal terms (16,17). 

The calculation of these terms can be subdivided in two steps. First, with help of a 
special subroutine we find the eigenvalues and eigenvectors of the shear matrix R (i. e. the 
quantities fm and the matrix O, see eq. (6)). The pressure tensor (17) is diagonal in the 
frame of eigenvectors of the matrix R. Therefore, we have to evaluate numerically only four 
quantities: three eigenvalues of the tensor (17) and the quantity (16). It turns out that in 
the case of polytropic stellar gas, the eigenvalues of the pressure tensor are proportional to 
the following integrals: 

f'2'K PIT g2 

where Cj are the direction cosines and 7 = 5/3 is the specific heat ratio. Quantity (16) is 
proportional to 

''^^ dd) r dOsinie) ^ ^, B2 

Integrals (Bl) and (B2) are evaluated numerically by a separate program and tabulated as 
functions of the ratios /(2)//(l) and /(3)//(l). For small and large values of the ratios we 
use analytical representations of these integrals in terms of series. 
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